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On this slide we show the three-dimensional unsteady form of the Navier-Stokes Equations These 
equations describe how the velocity, pressure, temperature, and density of a moving fluid are related. 
The equations were derived independently by G.G. Stokes, in England, and M Navier, in France, in the 
early 1800's. The equations are extensions of the Euler Equations and include the effects of viscosity on 
the flow. These equations are very complex, yet undergraduate engineering students are taught how to 
derive them in a process very similar to the derivation that we present on the conservation of momentum 
web page. 

The equations are a set of coupled differential equations and could, in theory, be solved for a given 
flow problem by using methods from calculus. But, in practice, these equations are too difficult to solve 
analytically. In the past, engineers made further approximations and simplifications to the equation set 
until they had a group of equations that they could solve. Recently, high speed computers have been 
used to solve approximations to the equations using a variety of techniques like finite difference, finite 
volume, finite element, and spectral methods. This area of study is called Computational Fluid 
Dynamics or CFD. 

The Navier-Stokes equations consists of a time-dependent continuity equation for conservation of mass, 
three time-dependent eOMervation of .momentum equations and a time-dependent conservation of 
energy equation. There are four independent variables in the problem, the x, y, and z spatial coordinates 
of some domain, and the time t There are six dependent variables; the pressure p, density r, and 
temperature T (which is contained in the energy equation through the total energy Et) and three 
components of the velocity vector; the u component is in the x direction, the v component is in the y 
direction, and the w component is in the z direction, All of the dependent variab/es are functions of ad 
four independent variables. The differential equations are therefore partial differential equations and not 
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the ordinary differential equations that you study in a beginning calculus class. You will notice that the 
differential symbol is different than the usual "d /dt M or "d /dx" that you see for ordinary differential 
equations. The symbol look like a backwards "6", and, because of font limitations, we will use the 
symbol "6 /6t" on this page to indicate partial derivatives. The symbol indicates that we are to hold all of 
the independent variables fixed, except the variable next to symbol, when computing a derivative. The 
set of equations are: 

Continuity: 6r/6t + 6(r * u)/6x + 6(r * v)/6y + 6(r * w)/6z = 0 

X - Momentum: 6(r * u)/6t + 6(r * u A 2)/6x + 6(r * u * v)/6y + 6(r * u * w)/6z == - 6p/6x 

+ 1/Re * { 6tauxx/6x + 6tauxy/6y + 6tauxz/6z} 

Y - Momentum: 6(r * v)/6t + 6(r * u * v)/6x + 6(r * v A 2)/6y + 6(r*v* w)/6z « - 6p/6y 

+ 1/Re * { 6tauxy/6x + 6tauyy/6y + 6tauyz/6z} 

Z - Momentum: 6(r * w)/6t + 6(r * u * w)/6x + 6(r * v * w)/6y + 6(r * w A 2)/6z = - 6p/6z 

+ I /Re * { 6tauxz/6x + 6tauyz/6y + 6tauzz/6z} 

Energy: 6Et/6l + 6(u * Et)/6x 4- 6(v * Et)/6y + 6(w * Et)/6z = - 6(r * u)/6x - 6(r * v)/6y - 6(r * w)/6z 

- 1/(Re*Pr) * { 6qx/6x + 6qy/6y + 6qz/6z} 

+ I /Re * {6(u * tauxx 4- v * tauxy + w * tauxz)/6x + 6(u * tauxy + v * tauyy + w * tauxz)/6y + 6(u * 

tauxz + v * tauyz + w * tauzz)/6z} 

where Re is the Reynolds number which is a similarity parameter that is the ratio of the scaling of the 
inertia of the flow to the viscous forces in the flow. The q variables are the heat flux components and Pr 
is the Prandtl number which is a similarity parameter that is the ratio of the viscous stresses to the 
thermal stresses. The tau variables are components of the stress tensor. A tensor is generated when you 
multiply two vectors in a certain way. Our velocity vector has three components; the stress tensor has 
nine components. Each component of the stress tensor is itself a second derivative of the velocity 
components. 

The terms on the left hand side of the momentum equations are called the convection terms of the 
equations. Convection is a physical process that occurs in a flow of gas in which some property is 
transported by the ordered motion of the flow. The terms on the right hand side of the momentum 
equations that are multiplied by the inverse Reynolds number are called the diffusion terms. Diffusion 
is a physical process that occurs in a flow of gas in which some property is transported by the random 
motion of the molecules of the gas. Diffusion is related to the stress tensor and to the viscosity of the 
gas. Turbulence, and the generation of boundary layers, are the result of diffusion in the flow. The Euler 
equations contain only the convection terms of the Navier-Stokes equations and can not, therefore, 
model boundary layers. There is a special simplification of the Navier-Stokes equations that describe 
boundary layer flows. 

Notice that all of the dependent variables appear in each equation. To solve a flow problem, you have to 
solve all fwe equations simultaneously; that is why we call this a coupled system of equations. There 
are actually some other equation that are required to solve this system. We only show five equations for 
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six unknowns. An equation of state relates the pressure, temperature, and density of the gas And we 
need to specify all of the terms of the stress tensor. In CFD the stress tensor terms are often 
approximated by a turbulence model. 
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Beginner's Guide Home Page 
NASA Glenn Learning Technologies Home Page 
http://www.grc.nasa.gov AVWW/K-1 2 



by Tom Benson 

Please send suggestions/corrections to: TIwttms JMemnn(ivjuisaj>o\> 
last Ufxlafeci Wed, Apr 14 03:55:22 PM EOT 2004 by jam Benson 
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